SHAPE MINIMIZATION OF DENDRITIC ATTENUATION 



ANTOINE HENROT AND YANNICK PRIVAT 

ABSTRACT. What is the optimal shape of a dendrite? Of course, optimality refers to 
some particular criterion. In this paper, we look at the case of a dendrite sealed at one end 
and connected at the other end to a soma. The electrical potential in the fiber follows the 
classical cable equations as established by W. Rail. We are interested in the shape of the 
dendrite which minimizes either the attenuation in time of the potential or the attenuation 
in space. In both cases, we prove that the cylindrical shape is optimal. 



1. Introduction 

1.1. Motivation. Is Nature always looking for optimum for living organisms? In par- 
ticular, are the organs designed to optimize some criterion? Complete answers to these 
questions are likely never to be discovered. Nevertheless, assuming that Nature proceeds 
in the most efficient way, can lead to a better understanding of the modeling of an organ 
and the underlying physical or chemical phenomena. This is this idea of inverse modeling 
that we had in mind when we began this work. Roughly speaking, it can be described by 
the following steps: 

i Let us consider a given organ of a living body. 

ii Write a mathematical model which describes the behavior of this organ. 

iii Imagine a (numerical) criterion that Nature would like to optimize for this organ. 

iv Determine the optimal shape for this criterion and this model. 

v Compare with the real shape(s). 

If the optimal and the real shapes coincide, we can guess that our model and our criterion 
are relevant. If they do not, we must admit that either our criterion or our model (or 
the initial guess that Nature looks for optimum) is probably wrong. We believe that it 
will often be the choice of the criterion which is not correct. A possible reason is the 
complexity of Nature. This complexity indicates that, in general, there is not a unique 
criterion to optimize but several ones (which could also be antagonists). The mathematical 
study (point 4 in the above procedure) becomes then much more difficult since one needs 
to use tools of multi-criteria optimization. 

In this paper, we want to follow the above procedure in the case of a dendrite. We consider 
a fiber which is sealed at its right end and connected to a soma at its left end. We use 
the classical cable equation to describe the electrical potential along the fiber. What are 
the criterions that we can consider here? Of course, we want the dendrite to propagate 
the best as possible the electrical signal. In other terms, the attenuation of the signal must 
be as small as possible. We are going to consider the two kinds of possible attenuation: 
attenuation in time or in space and we are looking for the shape of a dendrite which 
minimizes this attenuation. In both cases, the optimal shape that we get is a cylinder. Since 
it is very close to the real shape, at first sight, we can conclude that Nature is in accordance 
with mathematics for this problem and solves a shape optimization problem! For a general 
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reference in mathematical modeling in Neuroscience, we refer to the book of A. Scott, 
lfl5l . For a more exhaustive view and an introduction to the beauty of shapes (and, in 
particular, optimal shapes) in Nature, we refer to the classical books of S. Hildebrandt and 
A. Tromba, @ and A. Bejan, E). 



1.2. The mathematical model. Let us consider a fiber with a cylindrical symmetry, of 
length I and radius a(x) at point x. We denote by v(x, t) the difference from rest of the 
membrane potential at point x and time t. The equation satisfied by v(x, t) is similar to 
the classical cable equation as established by W. Rail during the sixties, cf lfl3l . lfT4l . 0. 
See also, lfl6l as the best motivation of the current study. We consider here the case of a 
fiber which is sealed at its right end and connected to a soma with surface area A s at its left 
end. Let us denote by R a the axial resistance (kilcm), C m is the membrane capacitance 
(/jF/cm 2 ), G m the fiber membrane conductance and G s the soma membrane conductance 
(mS/cm 2 ). We assume that the fiber is initially at rest and that it receives a transient current 
stimulus io at the left. The parabolic equation satisfied by v is then (see J4)): 



(1) 



f 2^^( fl2 i) =aVY+^P{C m f t +G m v) xe(0,£),t>0 

2^Mg2(0,i) = A s (C m f (0,t) + G s v(0,t)) - i (t) t > 

^(£,t)=0 t>0 

I v(x,0) =0 x e (0,f). 



It is convenient to represent v, solution of ([TJ, in terms of eigenfunctions as did S. Cox and 
J. Raol in 0): 

-f-oo 

(2) v(x,t) = J2M t )<t>n( x ) xe(0J),t>0 

n=Q 

where <fi n is the n-th eigenfunction associated to the eigenvalue [i n : 



-(a 2 <)' = ii n a^T+^$ n x e {0,1) 

^a 2 (O)0; i (O) + ( Mn + 7 )^(O) = O 


where 7 := 2R a (G m — G s ) is assumed to be non negative. We choose to normalize the 
eigenfunctions by 



(4) ||0„|| 2 :=A0 2 (O)+ / a(xWl + a' 2 (x)4> 2 n (x)dx = l 







where A = Of course the eigen-pair (^„,(/)„) strongly depends on the taper a(x) 
so we will often denote it by (/z„(a), cj> n (a)). The eigenvalue problem (O is not classical 
since the eigenvalue appears in the boundary condition, see section|2]for more precisions. 

Inserting the decomposition (f2]i in the equation (HJ gives an ordinary differential equation 
satisfied by tp n (t)- After resolution, we get: 



(5) v(x,t) = —— ^<t> n {0)<t> n {x)i *e- Xnt 



where X n := A1 " 2 t ? 2 " f ^° Gm > (see Lemma |2~TT ) and * denotes the convolution product of 



2R a C 

distributions 



1.3. The optimization problems. We need now to give a precise statement to the opti- 
mization problems presented in the introduction. In that purpose, we have to choose the 
functions we want to optimize and the class of functions a(x) in competition. Let us begin 
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with this last point. Since the fiber must not collapse, it is natural to assume a lower bound 
for the functions a(x), so we fix a positive constant ao and we impose: 

(6) Va;e[0,4 a(x) > a > . 

Now the minimal regularity needed for a is clearly, according to system ([T]) or (O that 
the derivative a' exists (at least almost everywhere) and is bounded, so we choose to work 
in the class of Lipschitz continuous functions which is often denoted by W 1,co (Q, £). At 
last, we also need to put a constraint on the "cost" for Nature to build a fiber. It seems 
reasonable to consider that this cost is proportional to the surface area of the fiber. This 
surface area is clearly given by 



(7) 



Surface area = 2n / a(x)\/\ + a' 2 (x)dx 
Jo 



so we can assume a bound, say S, on this surface area. To summarize, we consider the 
class of functions a(x) defined by: 



(8) A a0! s ■■= UieW 1 ' 00 ^,^, a(x)>a znd J a(x) y/l + a' 2 {x)dx < S 
Of course, we need to assume S > a^i in order that the class A aQ ,s be non trivial. 



As explained in the Introduction, an "ideal" dendrite should conduct in the best possible 
way the electrical information he is supposed to transmit. In other terms, the attenuation 
of the signal must be minimized. Since the potential v depends on the space and the time 
variable, we can consider both criterions. 

1.3.1. Attenuation in space. Let us introduce the transfer function T defined by : 

(9) T{a) := 7+°° ) Li 

Jo v(£,t)dt 

T(a) corresponds to the ratio of the mean values in time of the potential v taken at points 
x = and x = I. This ratio is always greater than one, see Remark [3~T1 and is a good 
indicator of the attenuation of the signal between the two extremities of the dendrite. So, 
it is realistic to look for a taper a{x) which yields a ratio as close to one as possible: 

(10) Find a in the class A ao ,s which minimizes T(a). 

1.3.2. Attenuation in time. According to expansion (0, the potential v(x,t) goes to 
when t — > +oo. More precisely, its asymptotic behavior is described by 

v(x,t)^-^—^ 1 {0) ( j> 1 (x)i o *e- Mt 

where Ai := ^^a" 2 ^ aGm > and ^ti(a) is the first eigenvalue of ©. Therefore, as it 
is classical in such parabolic problems, it seems natural to look for a function a(x) which 
minimizes the exponential rate of decay: 

(11) Find a € A a „,s which minimizes /ii (a) (the first eigenvalue of Q). 

The idea of minimizing eigenvalues of such Sturm-Liouville operators is a long story and 
goes back at least to M. Krein in ifTOl . see also Q for a review on such problems. 

1.3.3. The main result. We state in the following Theorem the main results of this paper 

Theorem 1.1. i The unique minimizer of the eigenvalue /ii(a) in the class A ao ,s I s 
the constant function a = ao. 
ii The unique minimizer of the criterion T{a) in the class A ao ,s is the constant func- 
tion a = <iq. 
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In other terms, for both criterions, the optimal shape of a dendrite sealed at one end and 
connected to a soma at the other end is the cylindrical one! 

In his thesis and in a foregoing paper, see ifTTIl . fl2l . the second author studies the case of a 
dendrite sealed at both ends. From a mathematical point of view, it changes the boundary 
conditions in ([TJ and © which become homogeneous Neumann boundary conditions at 
both extremities. The result he obtains is the same for the case of attenuation in space, but 
it is different for the attenuation in time. Actually, there is no existence of a minimizer 
for /xi(a) (as usual in the Neumann case /Uo(a) = and /ii(a) denotes the first non-zero 
eigenvalue). Moreover, he is able to exhibit minimizing sequences which would produce 
very strange dendrites! 

1 .4. Notation. The set of notation used in this paper is summarized in this section. 

H^ 1,oc (0, £) the set of Lipschitz continuous functions defined on the interval [0, £]. 
the class of functions defined by 

| a 6 W^iOj), a(x) > a , a(x)y/l + a ,2 {x)dx < sX. 
norm defined on the space of bounded functions L°°(0, £) by 

ll/lloo = SUp te[(W] \f(t)\. 

inner product defined for two continuous functions / and g by: 

< f,g > a := Af{0)g(0) + a{x)^/T+W(x)f{x)g{x)dx. 

||.|| a norm induced by <.,. > a . 

8 a completion of the space of continuous functions C([0, £]) for the norm || .|| a 

L 2 (0, £) the space of (classes of) functions which are square-integrable on (0, £). 

2 , . the Sobolev space of functions in L 2 (0, £) whose derivative (in the sense 

^ ' > ofdistributions)liesinL 2 (0,^). 

. dJ . Gateaux-derivative of a function J at point vq in direction h defined by: 



A 



■a ,S 



<-,->c 



2. Minimization of the first eigenvalue 

The eigenvalue problem (© is not completely classical due to the presence of the eigen- 
value //„ in the first boundary condition. As explained in (4), see also the works of J. 
Walter iflTI and J. Ercolano-M. Schechter 0, a good way to handle with such case con- 
sists in introducing the following inner product: 

< f,g> a = Af(0)g(0) + [ a(x) y/1 + a' 2 (x)f(x)g(x)dx, 
Jo 

its associated norm \ \.\\ a and the Hilbert space £ a defined as the completion of the space of 
continuous functions C([0,^]) for this norm. It is easy to see that £ a is a space satisfying 
i/ 1 (0,£) C £ a C L 2 (0,£) (both inclusions are strict). Moreover, the map i— ► 0(0) 
defines a linear continuous form on £ a . It is now classical spectral theory which allows to 
prove existence of a sequence of eigenvalues \i n and eigenfunctions 4> n orthogonal for the 
inner product < ., . > a . 

Let us now make an elementary observation on the sign of /xi (a): 
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Lemma 2.1. Let a be in the class A ao ,s> then the first eigenvalue /ii(a) o/([3| satisfies 

~2R a G m < -7 < m{a) < 0. 

Proof. Let a € A ao ,s an d i> € i? 1 (0,^), non identically zero. We denote by i?[a;i>], the 
Rayleigh quotient: 



R[a; v] 



a 2 (x)u' 2 (x)dx- A<yu 2 {0) 



a{x)xjl + a' 2 (x)u 2 {x)dx + Au 2 (0) 



The classical Poincare-Courant-Hilbert principle writes: 

(12) LLi(a) = inf R\a:v] 

veH 1 (o,l) 

Now, taking v = 1 in the above formula yields: 

A'y 

Mi (a) < Tt. < . 

a(x)y / l + a' 2 (x)dx + A 

The lower bound is easy by observing that \i\ (a) + 7 = R\a;<j>\ (a)] + 7 > 0. □ 
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Remark 2.1. Using the min-max formulae for the second eigenfunction, it is also possible 
to prove that the second eigenvalue satisfies /12 (a) > 0. The proof consists in studying the 
problem of calculus of variations min gL_r . J Q a 2 (x)v' 2 (x)dx on the class W a :— {v £ 

iJ^M) : v(0) /= Oand <u>i(a) > a = 0}. 

We can now prove the first part of Theorem 11.11 The eigenfunction 4>i (a) associated to 
Hi (a) realizes the minimum of the Rayleigh quotient. Hence, we have: 

rt 

a 2 {x)4> 1 {a)' 2 {x)Ax - A 7 0i(a) 2 (O) 

Mi (a) 



a{x)yjl + a' 2 (x)ct) 1 (ay(x)dx + ^i(a) (0) 
By Lemma l2TI the numerator of this quotient is negative. Moreover, we have: 

(13) / a 2 (x)Ma)' 2 {x)dx- A-y^a) 2 ^) > a 2 / ^{a)' 2 {x)dx - A 7 ^i(a) 2 (0) 
Jo Jo 

and 

1 1 

(14) - — . , , < 



/„ ay/ 1 + a' 2 Ma)(x)dx + Afafa) (P) /„ a o 0i(a)(a;)da; + A0i(a)^(O) 
(with a strict inequality in (fT3l l. ( TT4T > if a is not constant). Finally, writing that /Xi(ao) = 
inf i?[ao;u], we deduce from ( fT3b , (TT4T) and (TT2l) that Ui(ao) < ( a ) as soon as 

a /= a . 

3. Minimization of T(a) 

3.1. Introduction. This last section is devoted to the proof of the second claim of Theo- 
rem [TTT] i.e. that a = gq minimizes the criterion T(a). We recall that the criterion T(a) 
describes the attenuation in space and that it is defined by 

= J o + °°v(0,t)dt 

J +OC v(l,t)dt- 

The proof here is much more complicated than for \i\ (a). Let us now outline the different 
steps of the proof. 
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1st step: Using the Laplace Transform v(x,p) of v(x,t), we rewrite the criterion 
T(a) as the quotient v(0, 0)/v(£, 0). 

2nd step: We use the change of variable defined by y = -^-pj to transform the 
equation satisfied by v into a simpler differential equation. This allows us to con- 
sider a new unknown p(y) := a?(x)^J\ + a'(x) 2 instead of a and a new criterion 
T\ (p). The function p must lie in the set defined by: 



n 



a ,S 



:= j p G L°°(0, h) : al < p(y) and J ' p(y)dy < S 



3rd step: We solve the new optimization problem min T\ (p) first on the subclass of 
functions p G 7Z a „.s which satisfy p < M for some positive constant M, We prove 
that the minimizer has to be a bang-bang function. It means that it can only takes 
the values clq and M, Then, studying carefully the optimality conditions, we prove 
that the only minimizer is p = Oq. 

4th step: We conclude. 

3.2. Use of the Laplace Transform. The parabolic equation is not completely standard in 
the sense that it contains a dynamical boundary condition at x = 0. This kind of problem 
has been studied by different people, see e.g. J6), HI and the references therein. It can 
be proved that the solution v(x,t) belongs to L 2 (0, T, i? 1 (0, £)). Moreover, using the 
eigenvalue expansion (0, one can see that, in the case of an impulsion io = S (a Dirac 
measure at t — 0), the integrals Jq + °° |i>(0, t)\dt and Jq + °° \v(£, t)\dt are well defined. Let 
us introduce the Laplace Transform v(x, p) of the solution defined by 

r + co 

v(x,p) := / e- pt v(x,t)dt . 



o 



Thanks to the convergence of the integrals, the criterion T{a) can be rewritten as 

n ,, rpi s /o +0 °^(0^) d ^ lim p ^ ^(0,p) 

' W ' j+°°v{t,t)dt lim^o^rf' 

Now, transforming equation we see that the Laplace Transform v is the solution of the 
following o.d.e.: 
(16) 



1 d ( n 2dv 



a ~- = a 



VI + a' 2 {C m pv + G m v) (x,p) G (0,£) x (0;+oo) 



2R a dx V* dx 

|| ( ,p) = A s [C n pv(0,p) + G s v(0,p)] - 1 p G (0, +c^) 

I v(x,0) = x G (0,£). 

3.3. A change of variable. We are now going to use the following change of variable 
classical in ordinary differential equations, see e.g. O: 



V = 

The interval (0, £) becomes (0, l\) where 



At 



df 



« 2 W 

the function v is transformed into the function 

w(y,p) := v(x,p). 
and we consider a new unknown p defined by : 



(17) p(y) := a 3 (x)^l + a'(xy 
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Since a belongs to the class A ao ,s defined in (O, the new function p belongs to: 
(18) K ao , s := L G £°°(0,4) : a 3 < p(y) and jf ' p(y)dy < sj . 

Then, equation ( [Tol l becomes: 

2^0=P(CmP + G'm)w d/,p) G (0,4) x (0;+oo) 

^^(0,p)=A[CmP + G s mO,p)-l pG(0,+oc) 

%Vi,P) = pG(0,+oo). 

We let p going to in the equation ( fT9l (see AppendixlAlfor a mathematical justification) 
to get a function wo(y) := 0) which satisfies 



(20) 



m:^=pG m w yG(o,4) 



^^(0) = ^G sU ;o(0) 
^a(4) = 



Moreover, from (TT3T > the criterion T(a) becomes 
(21) T(a)=T l( p)- - o(0) 



The problems min{T(a), a G A aa .s} and mm{T\(p), p G 7Z ao .s} are n °t completely 
equivalent since a G *4 ao .s i— > p G 7t aa ,s is n °t a one-to-one correspondance. Neverthe- 
less, it is clear that lZ aa .s contains the whole image of A ao ,s by this map. So, if we find a 
minimizer of T\ in 7Z a0t s which belongs to the image of A aa ,s (this will be the case), we 
will solve our problem. 

Remark 3.1. Let us have a look to equation ( l20b . It is not possible that w' (Q) > (oth- 
erwise wo(0) would be positive, according to the first boundary condition, and then wq 
would remain positive and convex which would contradict the second boundary condition. 
In the same way, wo (0) cannot be negative, otherwise wo would remain negative and con- 
cave and this is impossible with the second boundary condition. So finally, one can see that 
w' Q (Q) < 0, wq(0) > and wq remains positive on the whole interval (0, At last, since 
w' is increasing and w' (£i) = 0, we see that wo is decreasing (and positive), therefore 
tu (4) < Wo(0) which proves that T (a) = Ti(p) < 1. 

3.4. Study of a new optimization problem. Using the different transformations intro- 
duced in the previous subsections, we see that we must now solve the optimization prob- 
lem: minTi(p) with p in the class lZ a(u s- Let us begin by solving this optimization prob- 
lem in the subclass 

Kl.s ■= \p G £°°(0,4) : 4 < p(y) < Mand ' p(y)dy < s\ 

where M is a fixed positive constant (M > Gq). We will let M — > +oo later. 

Theorem 3.1. The problem minTi(p), with p G s , has a solution p*. Moreover, 
every solution is a bang-bang function, i.e. a function which satisfies p* = afj or M 
almost everywhere. 



Proof. The fact that the optimizer is a bang-bang function often occurs in such control 
problems. Existence of a minimizer p* is easy, due to continuity of the criterion T\(p) for 
the weak-* convergence. Then, we write and analyze the optimality conditions thanks to 
the introduction of two adjoint problems. We are able to prove that the set {ag < p* < M} 
has zero measure. Let us now give the details. 
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Since TZ^f 5 is a bounded subset of L°°(0,4) it is compact for the weak-star conver- 
gence . So, to prove existence of a minimizer, we just need to prove that the criterion T\ 
is continuous for the weak-star convergence. Let (p n )neN be a sequence in Ti?£ s such 

that p n p and let us denote by w n and w the associated solutions of d40b . From the 
variational formulation of this problem 

(22) — / w' n (y)z'(y)dy + G m p n (y)w n (y)z(y)dy + 2AG s w n (0)z(0) = -±± 

LHa Jo JO 27T 

for every z <G H 1 (Q,£i), we first see (taking z = w n in (l22l i that the sequence w n is 
bounded in i7 1 (0,^i). So, it converges (up to a subsequence) weakly in H 1 (Q,£i) and 
strongly in L 2 (0,£i) to a function w^. Now, these convergence are sufficient to pass to 
the limit in 0221) . so we have proved that Woo = w and all the sequence converges since w 
is the only accumulation point. Existence of a minimizer p* in the class TZ^ s follows. 

We want now to write the optimality conditions. Let h £ W 1,oo (0, £\) be an admissible 
perturbation of the optimum. We will now denote by u>o the quantity: 



/dwo 
W0 := \~dp- (P U 

Classical variational analysis shows that u>o is the solution of the ordinary differential 
equation: 

' 2k^ L = G m (pw + hw ) ye (0,4) 
£^(0) = A s G s w (0) 

^(4) = o. 

Differentiating the criterion T\ in the direction h gives: 
(24) /fli A u)o(0)w (4) - w {0)w (h) 



Let us now introduce the two following adjoint problems. We consider the function q\ 
defined as the solution of the ordinary differential equation: 

2k^ = G mP ( qi (y)-y) ye (0,4) 
(25) { £^(0) = A s G sgi (0) 

and the function (72 solution of : 



2k^- = G mP (q2(y)-i) ye (0,4) 



2R a dy 



(26) \ ifrlf(O) = A s G s q 2 (0) 

^(4) = o 

(existence and uniqueness of q±, q^ follows from Lax-Milgram Theorem). Multiplying 
equation d23l l by </2 — 1, equation d26l i by wq and integrating both by parts yields: 

(27) «3o(0) = jf ' h(y)w (y)(q 2 (y) -l)dy. 

In the same way, multiplying equation d23l l by q\ — y, equation (l25l l by wq and integrating 
both by parts yields: 

(28) «5 (4) = w {0) + 2R a G m / h(y)w (y)(qi(y) - y) dy . 
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Therefore, d24l i together with d27l ) and (f28b gives 
(29) 

W' 7 = «F X %)wo(y) ( 1 (ft(v)-i)-«*(0)(9ifo)-y) 

where A = 2AG s R a . Let us denote by /, the function of one variable defined by: 

/: [0,4] — > R 



So, finally: 



dr \ r^ 1 

ft/ = J h(y)w {y)f(y)dy. 



We want now to prove that the optimum p* is a bang-bang function. It is a classical 
approach, see e.g. JS]. For that purpose, let us introduce the following sets: 



W) :={y€(0,£ 1 ):p*(y) = a }; 
l M (p*) := {y e (0,h) : p*(y) = M} ; 
2*(P*) : = ({V G (0,4) : «o < < M}. 



We write J^p*) = (J J y e (0,4) : a + -r < P *(y) < M - - I = (J X*, fc (p*). We 

fe=i ^ " J fe=i 

want to prove that T+^p*) has zero measure, for all integer k ^ 0. We argue by contra- 
diction. Let us suppose that \I+^{p*)\ ^ 0. Let y G l+^(p*)- We denote by {Gk, n )n>o, 
the sequence of subsets of ^ : 

G fc ,„ := B (</o, -) n2* >fc (p*) c Z* >fc (p*). 



Let us notice that |^| Gfc,„ = {yo}, an d let us choose /i = XG fc „■ Then, for t small 

n=0 

enough, perturbations p* + th et p* — th are admissible. Then: 

T^+tty-T^) = f h{y)M y )f{y)A y > o ^ / Mv)f(v)dV > 0. 

*V 1 JO JG fc ,„ 

In the same way: 

T^-th)-^*) = _ j h{y)My)f{y)dy > o ^ / WQ (y)f(y)dy < 0. 
f N° 1 JO JG t ,„ 

We can deduce that / wo(y)f(y)dy = 0. We divide by |Gfc, n | and we make n tending 

JG k , n 

to +oo. The Lebesgue density theorem shows that w(yo)f(yo) = 0, a.e. forj/o G Z*,k{p*)- 
This is clearly a contradiction, since wq and / are respectively solutions of the differential 
equations jRI^ = P*G m w Q and j^jjz = p*G m f. This proves that 
! i and then I+(p*) has also zero measure, what implies that p* equals af] or M almost 
everywhere. □ 

Now, we prove that among every bang-bang function, this is the constant function which 
yields the minimum of T\ . 

Lemma 3.1. The optimum ofT\ in the class s is the constant function p(y) = ag. 



Proof. Using notation of the proof of Theorem l3.ll the optimality conditions writes 

• On the setX (p*), we have f(y) > and h(y) > ; 

• On the set 1m{p*), we have f(y) < and h(y) < 0. 
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According to the differential equation (I26t and maximum principle, the function q2 — 1 is 
negative on [0, £ij. Let us write / like below: 

2R a G m f wo(h) - w (0) 



Wq(4)2 ™, , y A 

where g(y) := — — . The function g is two times derivable on (0,^i), and, using 
equations d25l l, d26l i we have: 

(30) Vy g (0,4), -7-n{y) = -2 —, >. — q— -^-(y)- 

dy 2 52(2/) -1 dy 

Therefore, on every interval where g| keeps his sign, we have: 

^-(y) = f 7 f 112 , with Cel. 
dy [92 (y) - 1] 

Now, is continuous on [0, £1], then the only possibility is: 

(3D V, e( 0, <1 ),|(„) = -||^>0. 



It follows that g is an increasing function on [0, £\] . Since 52 — 1 is negative and wq (0) > 0, 
we see that: /(y) > g(y) > t "^^~™ ( -°- ) . Then, according to the optimality 

conditions the function p* , local optimum for the criterion T\, is necessarily as follows: 

(32) p*{y) = { M 3 if f y< f 

y ' \ a 3 Q if y > 6 

with a transition point £1 which is possibly or t\ . Now, it remains to look for the better 



p(y) 



M 



Figure 1 . Possible profile of the optimum 



p* among all functions defined by (f32t . The only unknown is finally the transition point 
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£1. We can use the explicit expression of the solution wo of equation d40b for such a simple 
p*. We write wo 

i \ _ / a \ cosh(aj M y) + af sinh(w M y) on [0, £1] 
Wo[y > \ < cosh(w<,y) + sinh(w j/) on 



(33) 



3 — ™*ru»; 1 "-4 -v-ua; — i>± i ~±j 

~ four constants and ujq := ^/2R a G m aQ, ujm '■= \/2R a G m M. 
we get 



where af 1 , af, and af are iuui 
Thanks to boundary conditions, we get 

(34) to fa) = { a ™ cosh (™) + ^^'m 1 sinh (^My) on [0, £i] 

(cosh(woy) — tanh(wo^i) sinh(woy)) on [£1, £1 



Therefore, the criterion T\ is given by 



Ti(p*) 



w (0) af 



Finally, using the fact that wq is C 1 , we get, thanks to continuity of wq and at y = £1 

:h(w £i) cosh(w M £i) - — sinh(w A fCi) sinh(w £i) 
- tanh(cjo^i) sinh(cjoCi) coshfwM^i) sinh(a;M£i) cosh(wo£i) 



cosl 



(35) 
and then 

(36) 



cosh(wo^i) 
— sinh(wo^i) 



cosh(wo£i) cosh(wM$i) — sinhfwjvf^i) sinh(woCi) 

UJ M 

sinh(woCi) coshf^M^i) — sinhfwjvf^i) cosh(wo£i) 

UJ M 

Computing the derivative of the criterion with respect to the variable £1, we have 
(3?) d2^) = f<^) sinH ^ l)cosHLOo{ei _^ )] > , 

It follows that £1 has to be equal to 0, that means that the constant function Oq minimize 
the criterion T x on R*f „. □ 

3.5. Conclusion. The proof of Theorem 11.11 follows now easily. Since a§ is the unique 
minimizer of the criterion T\ in the class IZ^f s and 



s- 



M>a\ 



we get that a% is the (unique) minimizer of T\ in the class H ao ,S- Moreover, since o,q = 
a °\/ ^ (^f) ' * S c ^ ear tnat -^ 1 ( a o) = ^( a o) anc l men ' a o minimizes T on the class 

•Aa ,S- 



Appendix A. Limit of w p when p — > 

Let us now denote by w p the function w(.,p) for a given positive p. The equation ( fl9] i can 
be rewritten: 



1 d w p 



2R a 



-^f- = p (C m p + G m ) W p 



ye (0,4) 



(38) 



£^ (o) = A s [c mP + G s ] Wp (0) - 1 
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We recall that we want to prove that w p has a limit wq when p — ► and that wq is the 
solution of ( f20b - We can suppose that p G (0, 1]. Let us write the variational formulation 
of 01: 

Vz G i? (0, £i), a p (w p ,z) = L(z), 

with 

1 

otp(«,z) = — / u'(y)z'(2/) + (C m p + G m ) j o(2/)w( 2 /)z(y)d 2 / + A(C m p + G s )w(0)z(0) 

*tta Jo 

and 

Let us make z — w p in the above formulation. Since a p is uniformly (with respect to p) 
coercive: 



(39) 



iin (^-^ a o G m^j IWIhi(oA) < «p( w p. w p) = L K>) 



and L is a linear continuous form, we get from d39l that the sequence w p is bounded in 
iT^O^i). Therefore, there exists w* G H^iQ^iA such that w p w* in H^iQ^li) and 

w„ > w* in L 2 (0, ^i) (up to a subsequence). It follows that 

. tu p (0)— - *u;*(0), 

p->0 

p{y)w P {y)z{y)Ay — -> / p(y)w*(y)z(y)dy, 
Jo 

w' p {y)z(y)Ay -» / w*'(y)z'(y)dy. 

Therefore w* is the solution of the ordinary differential equation: 

2k^=pG m w J6(0,ii) 
^^(0) = ^G sU ;o(0)-l 
^(4) = 0. 
which means that = u>o, the desired result. 
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